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Abstract 

Arguments showing that exchange-only optimized effective potential (xOEP) methods, with 
finite basis sets, cannot in general yield the Hartree-Fock (HF) ground state energy, but a higher 
one, are given. While the orbital products of a complete basis arc linearly dependent, the HF 
ground state energy can only be obtained via a basis set xOEP scheme in the special case that 
all products of occupied and unoccupied orbitals emerging from the employed orbital basis set are 
linearly independent from each other. In this case, however, exchange potentials leading to the 
HF ground state energy exhibit unphysical oscillations and do not represent a Kohn-Sham (KS) 
exchange potential. These findings solve the seemingly paradoxical results of Staroverov, Scuseria 
and Davidson that certain finite basis set xOEP calculations lead to the HF ground state energy 
despite the fact that within a real space (or complete basis) representation the xOEP ground 
state energy is always higher than the HF energy. Moreover, whether or not the occupied and 
unoccupied orbital products are linearly independent, it is shown that basis set xOEP methods 
only represent exact exchange-only (EXX) KS methods, i.e., proper density- functional methods, if 
the orbital basis set and the auxiliary basis set representing the exchange potential are balanced to 
each other, i.e., if the orbital basis is comprehensive enough for a given auxiliary basis. Otherwise 
xOEP methods do not represent EXX KS methods and yield unphysical exchange potentials. The 
question whether a xOEP method properly represents a KS method with an exchange potential 
that is a functional derivative of the exchange energy is related to the problem of the definition 
of local multiplicative operators in finite basis representations and to the fact that the Hohenberg 
Kohn theorem does not apply in finite basis representations. Plane wave calculations for bulk 
silicon illustrate the findings of this work. 
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I. INTRODUCTION 



In a recent stimulating article with important implications for the use of finite basis 
sets, Staroverov, Scuseria, and Davidson^ presented an exchange-only optimized effective 
potential (xOEP) scheme that yields, for given finite Gaussian orbital basis sets, ground 
state energies that surprisingly equal exactly the ground state Hartree-Fock (HF) energies 
for these basis sets. Moreover, their xOEP scheme not only yields one unique but an infinite 
number of exchange potentials and each of the latter leads to the corresponding ground 
state HF energy if used as the exchange potential in the corresponding exchange-only KS 
Hamiltonian operator. On the other hand, it is known that in a complete basis set limit, 
which corresponds to a complete real space representation of all quantities, the xOEP method 
is identical^ to the exact exchange-only Kohn-Sham method and yields ground state energies 
that always lie above^ the corresponding ground state HF energy. Staroverov, Scuseria, and 
Davidson then state: "Our conclusions may appear paradoxical. For any finite basis set, no 
matter how large, there exist infinitely many xOEPs that deliver exactly the ground-state HF 
energy in that basis, however close it may be to the HF limit. Nonetheless, in the complete 
basis set limit, the xOEP is unique and E(xOEP) is above E(HF)". (Here E(xOEP) and 
E(HF) denote the xOEP and HF total energies, respectively, that are denoted E^^^^ and 
E^^ in this work.) Furthermore they state: "The non- uniqueness of OEPs in a finite basis 
set raises doubt about their usefulness in practical applications" 

We here first show, by different means including a constrained-search one, that the above 
statement of Staroverov, Scuseria, and Davidson, that it is always possible to construct 
optimized effective potentials that deliver exactly the ground state HF energy, holds if 
and only if the products of the orbital basis functions, or at least the products of the 
corresponding occupied and unoccupied HF orbitals from a given orbital basis set, form a 
linearly independent set. Otherwise, the xOEP scheme for finite orbital basis sets, in general, 
does not deliver exactly the ground state HF energy. Secondly, we show that the xOEP 
approach of Staroverov, Scuseria, and Davidson, does not really represent an exchange- 
only KS method and does not yield physically meaningfuU KS exchange potentials, even if 
the products of orbital basis functions are linearly independent. In order to get physically 
meaningfull KS exchange potentials via xOEP schemes, the latter have to be set up in 
a way that they represent KS methods, otherwise they are indeed of little usefulness in 
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practical applications. However, if xOEP schemes are set up properly then they are of 
great usefulness in practice as demonstrated, e.g., by numerically stable plane- wave xOEP 
procedures for solids^*^'^'^'^. 

II. RELATION OF xOEP AND HF ENERGIES WITHIN FINITE BASIS SET 
METHODS 

We start by briefly reconsidering the xOEP approach of Staroverov, Scuseria, and 
Davidson^. The relevant Hamiltonian operators are the HF Hamiltonian operator 

H^^ = -IV^ + vUr) + VH{r) + (1) 
and the exchange-only KS Hamiltonian operator 

= -l\/^ + v,^t{r)+VH{r)+iiA^) (2) 

Atomic units are used throughout. In Eqs. ([I]) and ([2]), fext(r) denotes the external potential, 
usually the electrostatic potential of the nuclei, vnir) is the Hartree potential, i.e., the 
Coulomb potential of the electron density, Vx{r) is the local multiplicative KS exchange 
potential, Vs{r) = vnir) + Vx{r) + fea;t(r) the effective KS potential, and v!^^ the nonlocal 
exchange operator with the kernel 

-^'M = ^ (3) 

in a real space representation. Here p(r, r') designates the first-order density matrix. In the 
HF-Hamiltonian operator of Eq. ([1]) the first-order density matrix occuring in the nonlocal 
exchange operator of Eq. equals the HF first order density matrix p^^(r, r') and the 
nonlocal exchange operator subsequently equals the HF exchange operator. For simplicity 
we consider closed shell systems with non-degenerate ground states. In this case orbitals, 
first-order density matrices, and basis functions can all be chosen to be real-valued. 

Next we introduce an orbital basis set {x^} of dimension A^. The representations of the 
HF- and exchange-only KS-Hamiltonian operators in this basis set are 

H^^ = T + V^^ + ^ + Ve.t (4) 
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and 



= T + Vn + Vx + Vext, (5) 

respectively. The matrices T, Yh, ^ext, ^x^ and are defined by the corresponding 
matrix elements T^^ = (x^l - IV^Ix,.), V^,^^ = {Xi,\vh\Xu), K^t.^i^ = {x^i\vext\Xl^) , V^^^ = 
{Xi^\v^^\Xu), and K,^^ = {Xi,\vx\Xu), respectively, and by = V// + + V^^f Because 
the orbital basis functions are real-valued all matrices are symmetric 

Now we expand the KS exchange potential in an auxiliary basis set {fk} of dimension 
M'^"^, i.e., 

^x(r) = ^ bkfkir) . (6) 
fe=i 

The auxiliary basis set, of course, shall be chosen such that its basis functions are linearly 
independent. The crucial question arising now is how many and what types of matrices V^; 
representing the KS density-functional exchange potential can be constructed for a given 
auxiliary basis set {fk}- This question was answered in Ref. 9j. Firstly we consider the 
case when the M = {1/2)N{N + 1) different products x^(r)xi,(r) of orbital basis functions 
are linearly independent. In this case, if M""^ = M and the auxiliary basis functions span 
the same space as the products of the orbital basis functions, then a symmetric matrix 
V^. can be constructed in a unique way by determining appropriate expansion coefficients 
bk for the exchange potential. The reason is that the determination of the M°"^ = M 
expansion coefficients bk for the construction of the M""^ = M different matrix elements of 
the symmetric matrix V^; leads to a linear system of equations 

Ab = W, (7) 

with 

^t,u,t = {Xf,Xv\fk) (8) 

for the coefficients bk of dimension M""^ = M that is nonsingular and thus has a unique 
solution^. In Eq. ([7]), A is a M x M""^ matrix that contains the overlap matrix elements 
{XtiXu\fk)- The first index of A, i.e., /zz/, is a superindex refering to products of orbital basis 
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functions while the second index, k, refers to auxihary basis functions. The vector b collects 
the expansion coefficients of Eq. ([6]) for the exchange potential and the right hand side W^;, 
a vector with superindices fiu, contains the M = N{N + l)/2 independent elements of an 
arbitrarily chosen matrix V^.. If we chose V^, to be equal to the matrix representation of an 
arbitrary nonlocal operator with respect to the orbital basis set then Eqs. (E]) and ([7]) define 
a local potential with the same matrix representation. This demonstrates that a distinction 
of local multiplicative and nonlocal operators is not clearly possible for orbital basis sets 
with linearly independent products of orbital basis functions. 

If M"""^ > M and the space spanned by the auxiliary functions contains the space spanned 
by the product of orbital functions then^ an infinite number of sets of coefficients bk lead to 
any given symmetric matrix V^.. The real space KS exchange potentials Vx{r) corresponding 
according to Eq. ([6]) to these sets of coefficients bk are all different but all represent local 
multiplicative potentials. Next we construct KS Hamiltonian operators ([2]) by adding these 
different KS exchange potentials to always the same external and Hartree potential. The 
resulting effective KS potentials in real space, i.e., the fs(r) are all different. Nevertheless 
the resulting basis set representations H^'^'^ of the corresponding KS Hamiltonian operators 
are all identical because the basis set representations V^. of the different exchange potentials 
Vx{r), by construction, are all identical. As a consequence the KS orbitals resulting from 
diagonalizing the KS Hamiltonian matrix H^'^"^ and subsequently also the resulting ground 
state electron densities are identical in all cases. We thus have a situation where different 
local multiplicative KS potentials fs(r) lead to the same ground state electron density. This 
seems to constitute a violation of the Hohenberg-Kohn theorem. Indeed it was shown in 
Ret. fl d.c..ed . Ret. S tha. .he Hohenbe..-Koh„ theo.en, doe= not hold to. 
finite orbital basis sets in its original formulation, i.e., that different local potentials, e.g., 
local potentials obtained by different linear combinations of auxiliary basis funtions, must 
lead to different KS determinants and thus different KS electron densities. We will come 
back to this point later on. Finally, if M"-"^^ < M then not all symmetric matrices \x can 
be constructed from a local KS exchange potential given by an expansion ([6]). 

In their xOEP approach Staroverov, Scuseria, and Davidson^ can expand the KS exchange 
potential in M""^ = M auxiliary basis functions and determine the coefficients such that 
the resulting matrix V^. exactly equals the HP exchange matrix V^^. If additionally the 
KS Hartree potential is set equal to the HP one then the resulting HP and KS Hamiltonian 
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operators are identical. Subsequently also the HF and KS orbitals, the ground state electron 
densities, and the ground state energies are identical. Because the HF and the KS electron 
densities turn out to be identical, the Coulomb potential of this density can equally well be 
considered as a HF or a KS Hartree potential. It follows immediately that the KS exchange 
potential constructed in this way is the xOEP exchange potential: The HF total energy 
is the lowest total energy any Slater determinant can yield. Thus if a local multiplicative 
KS potential leads to this total energy it is clearly the optimized effective potential defined 
as the potential that yields the lowest total energy achievable by any local multiplicative 
KS potential. The xOEP ground state energy resulting from this construction equals the 
corresponding HF energy. Moreover by enlarging the number of auxiliary basis functions, 
resulting in M"""^ > M, not only one optimized exchange potential leading to the HF energy 
but infinitely many can be constructed. 

Staroverov, Scuseria, and Davidson obtained the HF energy in their xOEP scheme even 
if the number of auxiliary functions only equaled the product M""" of occupied and virtual 
orbitals^. In this case a similarity transformation of the HF and the KS Hamiltonian matrices 
and their constituents was carried out in order to obtain representations of all matrices with 
respect to the HF orbitals. Then it is sufficient to chose the expansion coefficients of the 
KS exchange potential such that only the occupied- virtual block of the KS exchange matrix 
equals that of the HF exchange matrix. The resulting KS Hamiltonian matrix then may 
differ from the HF Hamiltonian matrix in the occupied-occupied and the virtual-virtual 
block but this merely leads to unitary transformations of the occupied and virtual orbitals 
among themselves and thus does not change the ground state energy or the electron density. 

Indeed it is straightforward to show that the occupied-virtual block of the exchange 
matrix equals that of the HF exchange matrix if the products of occupied and unoccupied 
orbitals are linearly independent. To that end we consider the xOEP equation determining 
the xOEP exchange potential^>^ 

occ. unocc. I I \ \ i \ occ. unocc. i , i Ml [ i \ 

4^ Y: Mr)Ur) = 4E E *,(r)0.(r) (9) 

I a I a 

In Eq. dHD (pi and (pa denote occupied and unoccupied KS orbitals, respectively, with 
eigenvalues Ei and Ea- Both sides of Eq. ^ are a linear combination of products 
0j(r)0a(r) of occupied and unoccupied KS orbitals with coefficients {(pal^xlfpi) I {^i — ^a) and 
{(t)a\v^^\(t)i) / {Ei — Ea), respectively. However, if the products 0i(r)0a(r) are linearly indepen- 

7 



dent then the two hnear combinations can only be identical if the coefficients multiplying 
the products are all identical. This, however, requires that (0a|fx|0j) = (^al'^i^^l^i), i-e., 
that the occupied-virtual block of the KS exchange matrix equals that of the corresponding 
exchange matrix of a nonlocal exchange operator of the form of the HF exchange operator. 
Replacement of the KS exchange matrix by the matrix of the nonlocal exchange operator 
thus again leads only to a unitary transformation of the occupied and virtual orbitals among 
themselves. Therefore the corresponding xOEP determinant can also be interpreted as HF 
determinant. 

Next we consider the crucial point what happens if the products of orbital basis functions 
X^(r)xj,(r) are linearly dependent. Then the rows of the matrix A of Eq. ([7]) are linearly 
dependent, thus the rank of the matrix A is lower than M, and as consequence Eq. ([7j), 
in general, has no solution. For an alternative argument, observe that for linear dependent 
products of orbital basis functions x^(r)xi,(r), there exists at least one linear combination 
of such products that equals zero 

= ^ af,^X,,{r)Xu{r) ■ (10) 

In Eq. (fTOl) the a^^iy denote the coefficients of that linear combination. The corresponding 
sum of matrix elements of V^; also equals zero, i.e., 

= ^ cif^u {Xf,\vx\Xu) = / drv^{r) ^ a^,^ Xi,{r)Xu{r) (11) 

for any choice of expansion coefficients ht in Eq. (I6l) because the product of any local function 
and thus of any KS exchange potential fx(i') with the sum ffTOj) equals zero. The products 
X^(r)Xi'(i"') foi' two different arguments r and r', on the other hand, are always linearly 
independent because the orbital basis set {Xa*} has to be linearly independent. Therefore 
the linear combination X^^i/ (^ij.vXii{^)Xu{^') can not be identical to zero for all values of 
the arguments r and r'. Then, however, also the integral of this linear combination with 
p^^{r, r')/|r — r'|, i.e., with the kernel of the nonlocal HF exchange operator, in general, is 
not equal to zero, i.e., in general 

O^Yl <'^^AX^.\€^\X.) ■ (12) 

Comparison of Eqs. (fTT!) and (fT2l) shows that, in general, the exchange matrices V^. and V^''' 
are different no matter how the expansion coefficients bk of the KS exchange potential, Eq. 
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([6]), are chosen. This demonstrates that, in general, neither the xOEP scheme of Staroverov, 
Scuseria, and Davidson^ nor any other leads to an xOEP Hamiltonian operator that equals 
the HF Hamiltonian operator when the orbital basis products are linearly dependent. If we 
consider the version of Staroverov, Scuseria, and Davidson's xOEP scheme that refers only 
to the occupied-virtual block of the xOEP and HF exchange matrices then by completely 
analogous arguments it follows that this scheme only works if the products of occupied and 
unoccupied HF orbitals are linearly independent. However, in general, if the products of 
occupied and unoccupied HF orbitals are linearly dependent then it is not possible to obtain 
the HF ground state energy via an xOEP scheme. 

A. Constrained-search analysis 

Before we discuss the question how products of basis functions can become linearly de- 
pendent for given orbital basis sets we elucidate the situation from a constrained-search^ 
point of view. We start with a constrained-search proof that the xOEP ground state en- 
ergy, E^^^^, must equal the HF ground state energy, E^^ , in their common finite orbital 
basis, when there is no linear dependence in the products of orbital basis functions. To 
accomplish this we appeal to the work of Harriman^^. He showed that only one first-order 
density matrix may yield any density generated by a given finite orbital basis whose basis 
products form a linearly independent set. This means that since an idempotent first-order 
density matrix uniquely fixes a corresponding single determinant, it follows that only one 
single determinant, constructed from a given finite orbital basis whose products are linearly 
independent, may yield a density that is constructed from this same basis. Consequently, 
with use of a common finite orbital basis set, the xOEP single determinant must equal the 
HF single determinant if there exists an effective KS potential Vg in Eq ([2]) such that the 
corresponding KS ground state density is the same as the Hartree-Fock density. That this Vg 
exists for the situation when the basis products are linearly independent, as discussed above, 
follows from Ref.jol and was shown in practice by Staroverov, Scuseria, and Davidson^. 

What happens when the products are not linearly independent? Due to the idempotency 
property of the first-order density matrix for a single determinant, a density generated from 
a given finite orbital basis could still generate a unique determinant if the basis products 
are linearly dependent, provided that this linear dependency is mild enough^, i.e., if the 
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products of occupied and unoccupied orbitals remain linearly independent. However, if 
the linear dependency of the basis product pairs is not sufficiently mild, then the situation 
changes dramatically in that more than one single determinant will yield the same density 
from a given finite basis set^^. In this case we do not have equality E^^^^ = E^^. Instead, 
we have inequality E^^^^ > E^^, which arises from the following contradiction. 

Assume that the xOEP determinant equals the HF determinant through 

respective optimizations in their common finite orbital basis set. Then it follows that their 
densities must be the same. But, from a constrained-search analysis^, the xOEP determi- 
nant ^^'^^P would yield this HF density and minimize, within this common basis, just the 
expectation value ($|T|$) of the kinetic energy, while the HF determinant yields this 
HF density and minimizes, within the common basis, the expectation value ($|T + Kel'^') 
of the kinetic energy plus the electron-electron repulsion energy. Here T denotes the many- 
electron kinetic energy operator, VJ=e the corresponding electron-electron repulsion operator, 
and $ Slater determinants that yield the HF density. (Equivalently, the xOEP determi- 
nant would yield the HF density and minimize {^\H — V^eel'^') while the HF determinant 
yields this HF density and of course minimizes Here H denotes the many-electron 

Hamiltonian operator.). Because the Slater determinants ^^'^^^ and <l>^^ minimize differ- 
ent expectation values, i.e., ($|T|$) and {^\T + V^el'^'), respectively, they are different, in 
general, and the inequality E^^^^ > E^^ applies for this common finite orbital basis case. 
However, there is only one possible determinant $ that yields the HF density from a given 
finite basis when the basis products are linearly independent or the extent of linear depen- 
dency is weak. In this case both minimizations yield this one Slater determinant simply 
because both minimization only run over one Slater determinant. Thus there is no contra- 
diction and the finite basis set conclusion of Staroverov, Scuseria, and Davidson follows in 
that the equality E^^^^ = E^^ applies. Hence we are now able to provide the resolution of 
the xOEP paradox- stated by Staroverov, Scuseria, and Davidson: For a finite basis set case, 
no matter how large the basis, E^^^^ equals E^^ provided that the basis products form a 
linearly independent set or the extent of linear dependence is sufficiently weak. However, in 
going from any starting finite basis set to the complete basis set limit, E^^^^ may become 
greater than E^^ somewhere along the way because as more and more basis orbitals are 
added to the finite basis set, the onset of sufficient linear dependency eventually occurs (see 
Appendix) . 
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We have provided an explanation for what might very well seem counterintuitive to the 
reader without knowledge of the analysis provided here. As one keeps adding more and 
more orbital basis functions, both E^^ and E^^^^ decrease and they continually remain 
equal to each other. Past a certain critical point in the addition of orbital basis functions, 
however, E^^ and E^^^^ may start to differ from each other and E^^ keeps decreasing 
while the behavior of E^^^^ depends on the chosen orbital basis set and it might actually 
be that E^^^^ rises! The latter behavior for example occurs if the exact HF orbitals as 
they correspond to a real space representation are themselves chosen as the basis set. If 
the basis set is restricted to the occupied HF orbitals, E^^ and E^^^^ are of course equal. 
If unoccupied HF orbitals are added to the basis set, E^^ remains unchanged at first. In 
contrast, beyond a certain point E^^^^ raises. The cause, of course, is the appearance 
of sufficient linear dependence at the critical point. (Ref. {l| does analyse certain linear 
dependency situations but the authors do not discuss the energy consequences for finite 
basis sets.) 



B. Creation of linear dependence 

Next we consider how products of orbitals basis functions become linearly dependent. As 
example we consider a plane wave basis set corresponding to a unit cell defined by the three 
linearly independent lattice vectors ai, a2, and as. The plane waves representing the orbital 
basis set {xg} then are given by 

Xc(r) ^ ^ (13) 

with 

G = £bi + mba + nbg (14) 

and 

£,n,m G Z and \G\<Gcut- (15) 

In Eq. dHl), bi, b2, bs denote three reciprocal lattice vectors defined by the conditions 
a£ • bm = 27r(5£m for £,m = 1,2,3. By Z the space of all integer numbers is denoted, Gcut 
denotes the cutoff that determines the size of the plane wave basis set, and fl stands for 
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the crystal volume. We have assumed before that basis functions are real-valued. This is 
not the case for plane waves. However, we can always obtain a real valued basis set by 
linear combining all pairs of plane waves with wave vectors G and — G to real- valued basis 
functions. This real valued basis set and the original complex-valued plane wave basis set are 
related by a unitary transformation that does not change any of the arguments of this paper. 
All arguments therefore are also valid for the complex- valued plane wave basis sets considered 
here and below. The number M of basis functions roughly equals (47r/3)Gj?„j(y/87r'^). The 
exact value of M depends on whether reciprocal lattice vectors G that lie in the immediate 
vicinity of the surface of the sphere with radius Gcut have lengths that are slightly larger or 
slightly smaller than Gcut- The relation 



shows that the products of plane waves of the orbital basis set are again plane waves of the 
same type with reciprocal lattice vectors G + G' that obey the relation G-l-G' < 2Gcut- Due 
to the latter relation the number of different products XgXc is about 8 times as large as the 
number of orbital basis functions, i.e., equals about 8A^. If > 15 then 8A^ < A^(A^+l)/2. 
In this case the number of different products of orbital basis functions is smaller than the 
number of products of orbital functions. Thus some products of orbital functions are equal 
and thus linearly dependent. For realistic systems the number of plane wave basis functions 
is much larger than 15. In a plane wave framework therefore xOEP and HF methods, in 
general, lead to different ground state energies with E^^^^ > E^^ . Results from plane 
wave xOEP and HF calculations for silicon discussed below illustrate this point. 

III. RELATION OF xOEP AND EXCHANGE-ONLY KS METHODS 

In this Section we show that the xOEP approach of Staroverov, Scuseria, and Davidson^ 
does not really correspond to an exact exchange KS method and does not yield a KS exchange 
potential, irrespective of whether or not the products of basis functions of the chosen orbital 
basis set are linearly independent. To this end we consider the xOEP or exact exchange 
(EXX) equation written in a form that slightly differs from that of Eq. iQ 



y-l gi(G+G' 



1 



XG+G'(r) 



(16) 



Vv 




(17) 
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The response function Xg in Eq. (fTTI) is given by 

X,(r,r') = aY^^Y. '^^^'^'^"^'^^^'^^'^^^'^^ • (18) 



Eq. f|T7|) can be derived in completely different ways. Firstly, following Refs. [lOj and 
one can consider the expression of the HP total energy and search for those orbitals that 
minimize this energy under the constraint that the orbitals are eigenstates of a Schrodinger 
equation with an Hamiltonian operator of the form 

The search for these orbitals is tantamount to searching the optimal effective potential 
^xOEP ^ therefore the name optimized effective potential method. The optimized effective 
potential v^^^^ can always be expressed as 

v^OEP^^^ ^ ^^^^(r) + vh{v) + v^iv) . (20) 
with the Hartree potential given as the Coulomb potential of the electron density generated 



by the orbitals. As shown in Refs. lOj and HI] the optimized effective potential v^^^^ is 



obtained if the exchange potential potential Vx of Eq. ( l20l) obeys the xOEP or EXX equation 

Alternatively the xOEP or EXX equation (fTTI) can be derived within an exact exchange- 
only KS framework. The Hamiltonian operator H^^^ of the exact exchange-only KS equa- 
tion is given by Eq. ([2]) with the effective KS potential 

Vs{r) = Vext{r) + vnir) + f^(r) . (21) 

The KS exchange potential in Eq. ( 1211) is defined as the functional derivative of the exchange 
energy 



= - E E I'd- I'd-' '^'^'^^'^t'-tf'^'^'^''^ (22) 



with respect to the electron density p, i.e, as 



PoUowing Ref. 20|j2l[ | we now exploit that according to the Hohenberg-Kohn theorem there 



exists a one-to-one mapping between effective potentials Vg and resulting electron densities 
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p. Therefore all quantities that are functionals of the electron density, here in particular the 
exchange energy, can be simultaneously considered as functionals of the effective potential 
Vs- Taking the functional derivative dE^/ 5vs{r) of the exchange energy with respect to the 
effective potential Vs in two different ways with the help of the chain rule yields 



The functional derivative 5p(r)/5fs(r') equals the response function f|T8l) and the right hand 
side of Eq. equals the right hand side of the xOEP or EXX equation (fT7I) . Furthermore 
the response function Xs is symmetric in its arguments for real valued orbitals. Therefore 
Eq. (p^ is identical to the OEP or EXX equation (fTTIl . This shows that the exchange 
potentials arising in the xOEP and the exact exchange-only KS schemes and subsequently 
the xOEP and the exact exchange-only KS schemes itself are identical. The xOEP or EXX 
equation can be derived in various ways within a KS framework^^. A crucial point, however, 
is that all derivations within a KS framework rely on real space representations in the sense 
that functional derivatives are taken within real space because the KS exchange potential 
is defined in real space as the functional derivative 6Ex/6p{r). Thus the above conclusion 
that the xOEP and the exact exchange-only KS schemes are equivalent holds only in real 
space, i.e., if all quantities are respresented in real space. Calculations, however, are usually 
carried out in basis sets and we will show next that in this case an xOEP and an exact 
exchange-only KS scheme, in general, are not equivalent. 

The xOEP or EXX equation f|T71) turns into the matrix equation 




(24) 



X 



t 



(25) 



with matrix and vector elements elements 




(26) 



a 




(27) 



and 




(28) 



a 
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if an auxiliary basis set {fk} is introduced to represent the response function, the exchange 
potential, and the right hand side of the EXX equation f|T71) . For simplicity we assume at 
this point that the auxiliary basis set is an orthonormal basis set. This is actually the case 
for plane wave basis sets but not for Gaussian basis sets. However, without changing the 
following arguments we can assume that we have orthonormalized any auxiliary Gaussian 
basis set. 

As long as the orbitals are represented in real space there is an infinite number of them 
and the summations over unoccupied orbitals in the response function (ITSl) and the right 
hand side of the xOEP or EXX equation ( |T7I) remains infinite and complete. For simplicity 
we assume that the considered electron system is either periodic and thus exhibits periodic 
boundary conditions or, in case of a finite system, is enclosed in a large but finite box with 
an infinite external potential outside the box. Then the number of orbitals is infinite but 
countable. As long as all orbitals are taken into account in the summation over unoccupied 
orbitals, the basis set representation of the exchange potential resulting from the basis set 
xOEP or EXX equation (|25l) becomes the more accurate the larger the auxiliary basis set 
and converges against the real space representation of the exchange potential and can be 
interpreted both as exact exchange-only KS or xOEP exchange potential. 

This changes dramatically if the orbitals are represented in a finite orbital basis set. 
Then, provided a reasonable orbital basis set is chosen, the occupied and the energetically 
low unoccupied orbitals are well represented. Most of the energetically higher unoccupied 
orbitals, however, are not represented at all simply because a finite orbital basis set can not 
give rise to an infinite number of unoccupied orbitals. Moreover, the energetically higher 
orbitals arising in a finite orbital basis set are quite poor representations of true unoccupied 
orbitals. Let us now concentrate on the representation of the response function. The inte- 
grals {4>a\fk\4>i) occuring in the matrix elements (126!) of the response function contain the 
three functions (pi, (pa, and fk- The occupied orbitals (pi have few nodes and thus are rela- 
tively smooth functions. The energetically low lying unoccupied orbitals still are relatively 
smooth, the higher ones however, with an increasing number of nodes and with increasing ki- 
netic energy become more and more rapidly oscillating. For smooth auxiliary basis functions 
/fc the integrals {(pa\fk\(pi) approach zero if they contain an energetically high unoccupied 
orbital (pa because the product of the smooth functions fk and (pi again is a smooth func- 
tion and the integral of this smooth product with a rapidly oscillatory unoccupied orbital 
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(f)a is zero due to the fact that any integral of a smooth with a rapidly oscillating function 
vanishes. This means that for matrix elements Xg^ke of the response function with two suffi- 
ciently smooth functions fk and fe the summation over unoccupied orbitals in Eq. (12^ can 
be restricted to unoccupied orbitals (pa below a certain energy depending on the smoothness 
of the involved auxiliary basis functions fk and xe- For sufficiently smooth functions fk and 
Xe the contributing unoccupied orbitals 0a thus are well represented in a finite orbital basis 
set. Therefore the matrix elements Xg^ke of the response functions are correct for indices 
k and i referring to sufficiently smooth auxiliary basis functions. For a more rapidly os- 
cillating auxiliary basis function nonvanishing matrix elements {4>a\fk\4'i) with energetically 
high unoccupied orbitals 0^ occur. The energetically high unoccupied orbitals (pa, however, 
are poorly described in the finite orbital basis set and moreover there are too few of them. 
Therefore the matrix elements Xg^ke of the response functions turn out to be wrong if at 
least one index refers to a more rapidly oscillating auxiliary basis functions. Indeed, if an 
auxilliary function fk oscillates much more rapidly than the energetically highest unoccupied 
orbitals (pa obtained for a given orbital basis set then all matrix elements {(pa\fk\<Pi) and thus 
all corresponding elements Xg^ke of the response function is erroneously zero. 

For a given auxiliary basis set, according to the above argument, a representation of 
the response function is correct only if the orbital basis set is balanced to the auxiliary 
basis set in the sense that it describes well unoccupied orbitals up to a sufficiently high 
energy. Otherwise an incorrect representation of the response function is obtained. The 
matrix representation of the response function like the response function itself is negative 
semidefinite. This is easily seen if a matrix element of the type {f\Xs\f) for an arbitrary 
function / is considered. Such a matrix element is obtained by summing up the contributions 
occuring in the summation over occupied and unoccupied orbitals in Eq. f|T8l) . Each single 
contribution and thus also the complete sum is nonpositive. Therefore an insufficient orbital 
basis set leading to too few energetically high unoccupied orbitals results in eigenvalues of 
the response matrix that have a too small magnitude. Solutions of the matrix equation 
( |25|) are given by the product of the inverse of the response matrix with the right hand 
side of the equation, i.e., by X^^t. If contains eigenvalues that are too small then 
the corresponding eigenvectors contribute with a too large magnitude to the solution of 
equation fl25l) . The eigenvectors with too small eigenvalues correspond to rapidly oscillatory 
functions. Therefore the resulting exchange potential exhibits rapidly oscillatory features. 
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This is exactly what is observed in the xOEP scheme of Staroverov, Scuseria, and Davidson^. 
If the response matrix even contains eigenvectors with eigenvalues that are erroneously zero 
then an infinite number of solutions arise of the matrix equation (1251) corresponding to an 
infinite number of exchange potentials, which yield, within the finite basis set, the same KS 
orbitals. 

Therefore if the auxilliary and the orbital basis sets are chosen unbalanced, e.g., if one 
chooses a too small orbital basis set for a given auxiliary basis set or a too large auxiliary 
basis set for a given orbital basis set, then the resulting response matrix is corrupted 
and no longer represents a proper representation of the response function in real space. 
In this case the xOEP scheme no longer represents an exact exchange KS scheme and 
the resulting exchange potential is unphysical and no longer represents the KS exchange 
potential. However, even in this case the xOEP scheme still is a proper optimized potential 
scheme in the sense that it yields a linear combination of auxiliary basis functions that results 
in the lowest total energy for this orbital basis set that can be obtained if the exchange 
potential shall be a linear combination of the auxiliary basis functions. While the resulting 
exchange potential is unphysical and does not resemble the KS exchange potential it obeys 
the above requirement of the xOEP scheme. The reason is that the arguments used for 
the xOEP derivation of the real space EXX or xOEP equation can also be used if orbital 
and auxiliary basis sets are introduced whereas no analogue to the DFT derivation exists 
anymore in this case. 



IV. EXAMPLES 



We now illustrate the arguments of the previous two Sections by specific examples. These 
examples also demonstrate that an auxiliary basis set that consists of all products of occupied 
and unoccupied orbitals is not balanced to the corresponding orbital basis set in the sense 
that a correct representation of the response function and a proper KS exchange potential 
can not be obtained for such an auxiliary basis set. Firstly a system of electrons in a box with 
periodic boundary conditions and an external potential equal to a constant is considered. 
The box shall be defined by corresponding unit cell vectors with i = 1,2,3. If the box, 
i.e., the unit cell vectors, become infinitely large then the system turns into an homogeneous 
electron gas. The KS eigenstates 0g of such a system are determined by symmetry and are 
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simple plane waves xg as they are given in Eq. (fT3ll . All plane waves with G vectors of a 
length smaller than some given constant Gp, i.e., with |G| < Gp, shall represent occupied 
KS orbitals, all plane waves with |G| > Gp represent unoccupied KS orbitals. The maximal 
length Gp of the vectors G of the occupied orbitals determines the Fermi level. For the 
orbital basis set as well as for the auxiliary basis set we chose plane waves, xg and /g, 
respectively, again as given in Eq. (fT3l) . Thus, for the considered system, arises the special 
case that each orbital basis function xg represents a KS orbital 0g- Obviously, the cutoff 
Gcut of the orbital basis set has to be chosen equal to or larger than Gp. 

The matrix representation of the response function in the considered case is diagonal 
with diagonal elements 

^ = 4 V ('^Gt'/Gl'^G'+G)(0G'+G|/G0G') 

''^^ ^ (IG'P- |G' + GP)/2 

^ V ^ (IG'P - |G' + GR ■ ^^^^ 

The auxiliary basis set shall be characterized by the cutoff radius Gcwf ' i-^-' auxiliary 
basis set shall consist of all plane waves /g with < |G| < G^^f. Note that the auxiliary 
function with G = that equals a constant function has to be excluded from the auxil- 
iary basis set because the xOEP or EXX equation in agreement with the basic formalism 
determines the exchange potential only up to an additive constant. A constant function 
would be an eigenfunction of the reponse function with zero eigenvalue. Now three cases 
can be distinguished: (i) If Gcwf — ^cut — Gp then the corresponding matrix elements 
-^s,GG of the reponse function are obtained with their correct value in a basis set calcula- 
tion with an orbital basis set characterized by the cutoff radius Gcut because all unoccupied 
orbitals 0g'+g occuring in the summation in Eq. [29] can be represented by the orbital ba- 
sis set. (ii) If Gcut — Gp < Gcuf — ^cut + Gp then for the matrix elements gg with 
Gcut — Gp < |G| < Gcut + Gp incorrect values are obtained because some of the unoccupied 
orbitals 0g'+g occuring for these matrix elements in the sum in Eq. [29] can not be repre- 
sented in the orbital basis set and therefore are not taken into account. Because all terms in 
the sum in Eq. [29] have the same sign the magnitudes of the resulting matrix elements gg 
are too small, (iii) If Gcut + Gp < G'^^f then the resulting X^ not only contains elements with 
a too small magnitude but additionally all matrix elements gg with Gcut + Gp < |G| are 
erroneously zero because all of the unoccupied orbitals 0g'+g occuring in the summation 
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in Eq. [29] can not be represented in the orbital basis set and therefore are not taken into 
account. 

If the auxihary basis set is chosen to be the space spanned by all products of occupied 
and unoccupied orbitals then it consists of all plane waves /g with < |G| < Gcut + Gp, 
i.e., Gcuf = Gcut + Gp- Thus the auxiliary basis set is chosen according to the above cases 
(ii). Therefore some of the resulting matrix elements X^ gg of the reponse function are 
incorrect. This demonstrates that an auxiliary basis set given by all products of occupied 
and unoccupied orbitals is not balanced with the corresponding orbital basis set. 

The considered system is special in that the right hand side of the xOEP or EXX matrix 
equation is zero due to the translational symmetry. Therefore also the resulting exchange 
potential is zero or more precisely equals an arbitrary constant. If the auxiliary basis set is 
chosen according to the above cases (i) and (ii) then a basis set calculation yields the correct 
exchange potential, i.e., zero or a constant. If the auxiliary basis set contains functions ac- 
cording to the above case (iii), however, then the xOEP or EXX matrix equation erroneously 
has an infinite number of solutions that equal a constant plus an arbitrary contribution of 
auxiliary basis functions with Gcut + Gp < |G|. The reason why the correct exchange poten- 
tial is obtained for an auxiliary basis set chosen according to the above case (ii) despite the 
fact that in this case the response function is already corrupted is that for the special system 
considered here the right hand side of the xOEP or EXX matrix equation is zero. Therefore 
any values for the diagonal elements Xs gg that differ from zero lead to the correct result. 
However, in general the right hand side of xOEP or EXX matrix equation is not equal to 
zero and then a response matrix with eigenvalues with erroneously too small magnitudes 
leads to a wrong exchange potential that exhibits too large contributions from those linear 
combinations of auxiliary basis functions that correspond to the too small eigenvalues of the 
response matrix. This is demonstrated in the following example. 

We consider plane wave xOEP calculations for bulk silicon carried out with the method 
of Ref. [41. The integrable singularity occuring in HP and xOEP exchange energies in 
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The lattice 



plane wave treatments of solids is taken into account according to Ref. 
constant was set to the experimental value of 5.4307A. The set of used k-points was cho- 
sen as a uniform 4x4x4 mesh covering the first Brillouin zone. In all calculations, all 
unoccupied orbitals resulting for a given orbital basis set were taken into account for the 
construction of the response function and the right hand side of the xOEP equation. EXX 
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pseudopotentials2^i2s with angular momenta / = 0, 1,2 and cutoff radii, in atomic units, of 
rf|j^Q = 1.8, r;?J^^ = 2.0, and r;?J^2 = 2.0 were employed. The pseudopotential with I = 1 was 
chosen as local pseudopotential. 

Figs.[T]and[2]display xOEP exchange potentials along the silcion-silicon bond axis, i.e., the 
unit cell's diagonal, for auxiliary basis set cutoffs E^^^ of 5.0 and 10.0 a.u. (6*""^ 3.2 and 
4.5 a.u.), and for various different orbital basis set cutoffs Ecut- Note that in figures and tables 
instead of the cutoffs Gcut and Gcwf ^^at refer to the length of the reciprocal lattice vectors 
of the plane waves the corresponding energy cutoffs Ecut = ^G^^^ and E^^^ = )^ 
displayed. Fig. [1] shows that the combination of an auxiliary basis set with cutoff E^^t = 5.0 
i^cut — 3.2) with a orbital basis set with cutoff Ecut = 1-25 {Gcut = 1-6) leads to a 
highly oscillating unphysical exchange potential. The cutoff of the auxiliary basis set in the 
considered case is about twice as large as the cutoff of the orbital basis. This means that 
the space spanned by the auxiliary basis is the same as that of all products of occupied and 
unoccupied orbitals. In this case the matrix representing the response function is corrupted 
and the resulting exchange potential turns out to be unphysical. With increasing cutoff 
Ecut of the orbital basis set the xOEP exchange potentials converge towards the physical KS 
exchange potential, more precisely towards the representation of the physical KS exchange 
potential in an auxiliary basis set with cutoff -E^uf — 5-0 {Gcut — 3-2). If the cutoff Ecut of 
the orbital basis set is about 1.5 times as large as the cutoff of the auxiliary basis set E^^^, 
i.e., equals 7.5 {Gcut = 3.9), then the exchange potential is converged. A further increase of 
Ecut to Ecut = 10.0 {Gcut = 4.5) leads to an exchange potential that is indistinguishable from 
that for Ecut = 7.5 {Gcut = 3.9) on the scale of Fig. [H Fig. [2]gives an analogous picture for a 
cutoff of the auxiliary basis set of Ecut = 10-0 (^c«f = 4.5). Again, if the space spanned by 
the auxiliary basis set equals that of the product of occupied and unoccupied orbitals, curve 
for Ecut = 5.0 {Gcut = 3.2), an highly oscillating unphysical exchange potential is obtained. 
If Ecut ~ 1-5-E'cuf then the exchange potential is converged towards the representation of the 
physical KS exchange potential in an auxiliary basis set with cutoff Ecut = 10-0 (^cuf = 4.5). 

This demonstrates the point that the xOEP scheme only represents a KS scheme if the 
orbital basis set is balanced to the auxiliary basis set. In the case of a plane wave basis set 
this requires the energy cutoff Ecut of the orbital basis set to be about 1.5 times larger than 
the energy cutoff E^^^ of the auxiliary basis set. 

Table [T] lists for a number of orbitals basis set cutoffs Ecut exchange and ground state 
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energies for series of auxiliary basis set cutoffs E^^f. Table [T] shows that the ground state 
energies for a given Ecut always decrease with increasing E^^f even if the values of E^^f 
is that large that the resulting exchange potential is unphysical. This demonstrates that 
the xOEP scheme remains well-defined even if unbalanced basis sets are used. In this case, 
however, the xOEP scheme no longer represents a KS method and the resulting exchange 
potential is unphysical and does not represent the KS exchange potential. Table [T] also 
lists the differences of the xOEP and HF ground state energies and shows that the xOEP 
energy does not converge to the HF energy. In the combinations E^ut = 2.5/-E'^^f = 10.0, 
Ecut = ^-O/E^ut = 20.0, and Ecut = 7-5/E^^t = 29.9 the space spanned by the auxihary 
basis set roughly equals that of the product of occupied and unoccupied orbitals. The ground 
state xOEP energies in these cases is de facto the lowest that can be achieved by the xOEP 
method for the given orbital basis set. The fact that this energy is higher than the HF 
total energy shows that the xOEP energy does not reach the HF ground state energy if the 
products of occupied and unoccupied orbitals become linearly dependent as it is usually the 
case in plane wave calculations and as it is the case in the presented calculations. 



V. SUMMARY 



We have given arguments leading to the conclusion that exchange-only optimized po- 
tential (xOEP) methods, with finite basis sets, cannot in general yield the Hartree-Fock 
(HF) ground state energy, but a ground state energy that is higher. This holds true even 
if the exchange potential that is optimized in xOEP schemes is expanded in an arbitrarily 
large auxiliary basis set. The HF ground state energy can only be obtained via an xOEP 
scheme in the special case that all products of occupied and unoccupied orbitals emerging 
for the orbital basis set are linearly independent from each other. In this case, however, 
exchange potentials leading to the HF ground state energy exhibit unphysical oscillations 
and do not represent Kohn-Sham (KS) exchange potentials. These findings solve the seem- 
ingly paradoxical results of Staroverov, Scuseria and Davidson^ that certain finite basis set 
xOEP calculations lead to the HF ground state energy despite the fact that it was shown^ 
that within a real space representation (complete basis set) the xOEP ground state energy 
is always higher than the HF energy. A key point is that the orbital products of a complete 
basis are linearly dependent. 
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Moreover, whether or not the products of occupied and unoccupied orbitals are hnearly 
independent, we have shown that basis set xOEP methods only represent exchange-only 
(EXX) KS methods, i.e., proper density-functional methods, if the orbital basis set and 
the auxiliary basis set representing the exchange potential are balanced to each other, i.e., 
if the orbital basis set is comprehensive enough for a given auxiliary basis set. Otherwise 
xOEP schemes do not represent EXX KS methods. We have found that auxiliary basis 
sets that consist of all products of occupied and unoccupied orbitals are not balanced to 
the corresponding orbital basis set. The xOEP method, even in cases of unbalanced orbital 
and auxiliary basis sets, works properly in the sense that it determines among all exchange 
potentials that can be represented by the auxiliary basis set the one that yields the lowest 
ground state energy. However, in these cases the resulting exchange potential is unphysical 
and does not represent a KS exchange potential. Therefore the xOEP method is of little 
practical use in those cases for which it does not represent a EXX KS method. Remember 
that, at present, the main reason to carry out xOEP methods in most cases is to obtain a 
qualitatively correct KS one-particle spectrum, either for the purposes of interpretation or 
as input for other approaches like time-dependent density-functional methods. However, the 
unphysical oscillations of the exchange-potential of xOEP schemes with unbalanced basis 
sets affect the unoccupied orbitals and eigenvalues. Another reason to carry out xOEP 
methods that represent EXX KS methods is that the latter may be combined with new, 
possibly orbital-dependent, correlation functionals to arrive at a new generation of density- 
functional methods. Also in this case it is important that the xOEP methods represents 
proper KS methods. 

A balancing of auxiliary and orbital basis sets is straightforward for plane wave basis 
sets. In this case xOEP schemes are proper EXX KS methods if the energy cutoff for the 
orbital basis set set is about 1.5 times as large as that of the auxiliary basis set. This as well 
as other results of this work were illustrated with plane wave calculations for bulk silicon. 
For Gaussian basis sets on the other hand, a proper generally applicable and reasonably 
simple balancing scheme of orbital and auxiliary basis sets is so far not available despite 
much effort s2^i2£i2^i2^'^'^. Therefore effective exact exchange-only methods like the KLI— , 
the 'localized Hartree-Fock'— , the equivalent 'common energy denominator approximation' 
method^^, or the closely related very recent method of Ref. IsJ, are in use as numerically 
stable alternatives that yield results very cose to those of full EXX KS methods. 
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TABLE I: xOEP exchange and ground state energies and , respectively, and the 

difference between HF and xOEP ground state energy for various combinations of and orbital 
auxiliary basis sets, characterized by energy cutoffs Ecut and -Epuf > respectively. {N and 
denote the corresponding number of basis functions.) All quantities are given in a.u. 
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VI. APPENDIX: LINEAR DEPENDENCE OF PRODUCTS OF BASIS FUNC- 
TIONS OF A COMPLETE BASIS 

Let {(f)k{x)} be a complete set of functions of a complex valued variable x such that any 
arbitrary square integrable function can be written as a linear combination of the functions 
in the complete set. We show that the set {(j)k{x)(j)i{x)} is linearly dependent. 

Using our complete sets, an arbitrary function f{x,y) of two complex valued variables x 
and y may be expanded in terms of {^^(a;)} and {(t>i{y)} 

oo oo 

/(a;, y) = ^Yl ^^^^ My)Mx) (30) 

fc=l £=1 

Set y = X to get: 

oo oo 

f{x, a;) = ^ ^ bk,i 0^(x)0fc(x) (31) 

k=l 1=1 

Now choose a function f{x,x) and a 0„(a;) out of the set {(f)k{x)} such that (i) 
lim^^oo = ^^(^ (ii) least one bk/ 7^ when £ ^ n and k ^ n. Since ^^^0^ is 

just a function of x, we may expand it in term of the {0fc(x)}: 

^-^ = f2dm^^^^) (32) 

Solving for /(x, x), 

00 00 

/(X, X) = (l)n{x) ^ (im0m(x) = ^ rf„0^(x)0„(x) (33) 
m=l m=l 

and equating Eq. (13T1) with Eq. (!33l) . we get 

oo oo oo 

(im0„(x)(/)„(x) = ^ ^ bk/(pk{x)(l)e{x) (34) 

m=l fc=l £=1 

or by setting k = m, 

oo oo oo 

^ c?^0„(x)0„(x) - &m,£0m(a;)0Ka;) = (35) 

m=l m=l i=l 
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or 



oo 



oo 



oo 



•n 




)(j)n{x)(t)j{x) - ^ ^ hm,i(l)m{x)(l)t{x) = . (36) 



e=i m=l 



Eq. is a linear combination of a subset of {(f)k{x)(f)i{x)} broken up into disjoint com- 
ponents and equated to zero. If a subset of a set is linearly dependent, then the set must 
also be linearly dependent. We show such a case by contradiction: According to Eq. (l36l) . 
for the subset {(j)k{x)(j)i{x)}(^ (appearing in the equation) to be linearly independent, three 
conditions must be met: 



2. dj - bj^n - bn,j = (V j G N with j ^ n) 

3. bm,t = (V m, ^ G N with m ^ n and i ^ n) 

But according to our condition on /(x, x) there is at least one bm,t 7^ with m ^ n and 
£ 7^ n, which is a contradiction to number three of our linear independence criteria. Therefore 
{0fc(x)0i(a;)}c must be linearly dependent by contradiction, and therefore {(t)k{,x)(l)i{x)} for 
all A;, / G N is linearly dependent because {0fc(x)0;(x)}c is linearly dependent. 

One may take the result one step further to show with an induction argument that for 
any complete set such as {(j)k{x)} the set defined by {Hill </'pi(^)l ^^hPi ^ N} is complete 
and hnearly dependent. 



1. d, 
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[111] direction 

FIG. 1: xOEP exchange potential along the silcion-silicon bond axis, i.e., the unit cell's diagonal, 
for an auxiliary basis set cutoff = 5.0 a.u. and different orbital basis set cutoffs Ecut- The 

upper and lower panels differ in the energy scale. The curve for Ecut = 1-25 a.u. is only displayed 
in the upper panel. 
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[111] direction 



FIG. 2: xOEP exchange potential along the silcion-silicon bond axis, i.e., the unit cell's diagonal, 
for an auxiliary basis set cutoff = 10.0 a. u. and different orbital basis set cutoffs Ecut- The 
upper and lower panels differ in the energy scale. The curve for Ecut = 5.0 a.u. is only displayed 
in the upper panel. 
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